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Abstract 

In this paper, the global well-posedness and stability of classical solutions to the multi- 
dimensional hydrodynamic model for semiconductors on the framework of Besov space are 
considered. We weaken the regularity requirement of the initial data, and improve some 
known results in Sobolev space. The local existence of classical solutions to the Cauchy 
problem is obtained by the regularized means and compactness argument. Using the high- 
and low- frequency decomposition method, we prove the global exponential stability of classi- 
cal solutions (close to equilibrium) . Furthermore, it is also shown that the vorticity decays to 
zero exponentially in the 2D and 3D space. The main analytic tools are the Littlewood-Paley 
decomposition and Bony's para-product formula.. 
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1 Introduction and Main Results 

In the modern semiconductors industry, the numerical simulation of device modeling has become very 
important. With the fast development of miniaturization devices, the traditional drift-diffusion model is 
no more valid, especially in submicron devices or in the occurrence of high field phenomena. Some kinetic 
models such as Boltzmann equation which describe the evolution of the distribution function f(t, x, v) of 
charged particles in the phase space are more accurate, but they need much computing power in practical 
application. Therefore, the hydrodynamic model which represents a reasonable compromise between the 
physical accuracy and the reduction of computational cost has recently received increasing attention in 
statistical physics and applied mathematics. By applying a moment method and appropriate closure 
conditions, it can been derived from the semiclassical Boltzmann equation coupled with the electric 
potential through a Poisson equation (see, e.g., Refs. [H]). In this paper, we are interested in the 
simplified hydrodynamic model where the energy equation is replaced by a pressure-density relation, 
which was first analyzed by Degond and Markowich [7]. The main objective is to study the global well- 
posedness and stability of classical solutions to the (unipolar) hydrodynamic model on the framework of 
Besov space. After an appropriate scaling, it can be written as 

n t + V • (rau) = 

u t + (u • V)u + ±Vp(n) = V$ - h (1.1) 
A$ = n — n, $ — > as \x\ — > +oo 
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for (£, x) G [0, +00) x Mr, N ^ 2, where n, u = (it 1 , u 2 , • • ■, u N ) T ( T represents transpose) and $ denote 
the electron density, the electron velocity and the electrostatic potential respectively. V is the gradient 
operator and A is Laplacian operator. The constant t is the momentum relaxation time of electrons. 
The pressure p satisfies the usual 7-law: 

p = p(n)=An 1 , (1.2) 

where the case 7 > 1 corresponds to the isentropic gas and 7=1 corresponds to the isothermal gas, A 
is a positive constant. The constant n > stands for the density of positively charged background ions. 
The system is supplemented with the initial data 

(n,u)(x,0) = (n ,u )(x), x G R N . (1.3) 

There are many contributions in mathematical analysis for (1.1)-(1.2), like well-posedness of steady 
state solutions, global existence of classical or entropy weak solutions, large time behavior of classical 
solutions, relaxation limit problems and so on, we may refer to Refs. [TJ [7l [lOl [T4J [18] and the 
references therein. For the evolutionary system (1.1)-(1.2), Luo, Natalini and Xin [M] first established 
the global exponential stability of small classical solutions to the Cauchy problem. Subsequently, there are 
many well-posedness and stability results for the Cauchy problem and the initial boundary value problem 
in one dimensional or multidimensional space. Physically, it is more important and more interesting to 
study (1.1)-(1.2) in the multidimensional case. However, up to now, only partial results are available. 
Hsiao and Wang et al. (TOj [18] studied the spherically symmetrical solutions to (1.1)- (1.2) with 7=1, 
O = {x G R N \ < R t < |sc| < R 2 < +00} and 7 > 1, fi = {x G \x\ ^ R 1 > 0} respectively. Guo 
investigated the irrotational Euler-Poisson equation (N — 3) without relaxation and constructed the 
global smooth irrotational solutions based on the Klein-Gordon effect, which decay to the equilibrium 
state uniformly as (1 + t)~ p (l < p < 3/2). Recently Hsiao, Markowich and Wang [5] dealt with the 
multidimensional unbounded domain problem (N — 2,3) without any geometrical assumptions. The 
main ingredient was to arrive at an a-priori estimate in terms of the classical energy argument. Later, 
Ahp] discussed the extended thermodynamic model (N 2) and reached the global existence, uniqueness 
and stability of classical solutions. The key a-priori estimate was obtained by the positive definiteness of 
some Liapunov functions. The above results are established on the framework of Sobolev space H (M ) 
and the regularity index is required to be high (£ > 1 + y > I G Z) when one deals with them by classical 
analysis methods. To lower the regularity, using the Littlewood-Paley decomposition and Bony's para- 
product formula, we prove the following well-posedness results for the system (1.1)-(1.3) in the critical 
nonhomogcncous Besov space 

B li(^ N ) = 1 + T )• First of all > 

we give a local existence and uniqueness 
theorem of classical solutions to (1.1)-(1.3) away from the vacuum. 

Theorem 1.1. Let n > be a constant reference density. Suppose that no — n, Uq and eo G B 1 ^ i(^ ) 
with ng > 0, then there exist a time To > and a unique solution (n, u, V$) of the system (l.l)-(l.S) 
such that 

(n, u, V$) G £^([0, T ] x R N ) with n>0 for all t G [0, T ] 

and 

(n - n, u, V$) G C([0, T Q ],B° tl {R N )) n ^([0, T ], B^ X (R N )), 
where eo := VA~ 1 (no — n). 
Remark 1.1. The symbol VA^ 1 means 

VA"V= f V x G(x - y)f(y)dy, 
Jr n 

where G(x, y) is a solution to A x G(x, y) = 5(x — y) with x, y G Mr. 

Remark 1.2. (1) When one estimates the nonlinear pressure term by the spectral localization method, 
there appear many commutators, which make computation very tedious. Thanks to the ideas in Ref. [16j . 
we introduce a function transform (sound speed) to reduce the system (1.1)-(1.3) to a symmetric hyper- 
bolic system (3.1)-(3.2) where the nonlinear term becomes linear and bi-linear terms. But, under the 
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transform there is a new function h(m) in Poisson equation, which is well-defined and smooth on the 
domain {m|^— m + ip > 0}. Theorem 1.1 follows from Proposition 4.1 and Remark 3.1. We extend the 
local well-posedness result in H e (M. N ) (see Ref. [1] or [9]) to that in Besov space m virtue of 

the regularized means and compactness argument. 

(2) Although our local existence result is proved via a symmetric hyperbolic system, there are some 
especial contents from Poisson equation to be dealt with. Therefore, we can not use the result of Iftimic 
[IT] directly. 

Under a smallness assumption, we establish the global existence, uniqueness and exponential decay 
of classical solutions to (1.1)-(1.3). 

Theorem 1.2. Let n > be a constant reference density. Suppose that uq — n, Uq and e$ G Sf lC^)- 
There exists a positive constant So depending only on A, 7, t and n such that if 

|| (n - n, it , eo)|| B j i(R iv) < S , 

then there exists a unique global solution (n, u, V$) of the system (1.1)- (1.3) satisfying 

(n,u,W) S ^([0,00) x R N ) 

and 

(n - ft, u, V$) e C([0, 00), ^ N )) n C\[0, 00), B^(R N )). 
Moreover, we have the decay estimate 

\\(n- n, ■u,V$)(-,i)|| B . i(R iv ) + || (?i t , it t , V$ t )(-, t)|| B j-i (RJV) 

< C 1| (no - n, uq, e )|| B j i(R iv) exp(-^ t), t > 0, 
where the positive constants /io and Co depend only on A, 7, r and n. 

Remark 1.3. From the proof of Theorem 1.2, one can see that the decay rate fiQ satisfies fio — Kt if 
< t ^ 1 and /Jo = K/t if t > 1, where if is a positive constant independent of t. 

Remark 1.4. Theorem 1.2 follows from Proposition 5.1 and Remark 3.1. Proposition 5.1 relies on a crucial 
a-priori estimate in Proposition 5.2 and the standard boot-strap argument. In the proof of a-priori 
estimate, we use high- and low- frequency decomposition method. Different from the classical energy 
argument in Ref. [9], our method shows that Poisson equation plays a key role in the low frequency 
estimates (see (5.3), (5.4) and (5.10)). Such fact leads to the global exponential stability of classical 
solutions to (3.1)-(3.2). 

Based on Theorem 1.2, we can characterize the exponential decay of the vorticity u) = V x u in Besov 
space BS^QR*). 

Theorem 1.3. (N=2,3) Let (n, u, V$) be the solution in Theorem 1.2. If 

||(n - n, uq, eo^lsg-^RW) < 
then the vorticity w — V x u decays exponentially in B^ 1 (R N )(u;(x, 0) = V x uq): 

IK'» *)lls-- 1 (M^) < ll W (-'°)lls-7 1 (R«) ex p(- 7)' i >°> 

where the positive constant 6' — min{(5o, 2Cr 1 CoT } depends only on A,j,t and n (C5 a constant given in 
(5.22)). 

The paper is arranged as follows. In Section 2, we present some definitions and basic facts on 
the Littlewood-Paley decomposition and Bony's para-product formula. In Section 3, we reformulate the 
system (1.1)-(1.3) in order to arrive at the effective a-priori estimates by the spectral localization method. 
In Section 4, we obtain the local existence and uniqueness of classical solutions to (3.1)-(3.2) with general 
initial data. In the last section, we deduce a crucial a-priori estimate under a smallness assumption, 
which is used to achieve the proof of global existence. Furthermore, it is also shown that the vorticity 
decays to zero exponentially in the 2D and 3D space. 

Throughout this paper, the symbol C denotes a harmless constant and all functional spaces are 
considered in M. N , so we may omit the space dependence for simplicity. Moreover, the integration J RN fdx 
is labeled as / / without any ambiguity. 
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2 Littlewood-Paley Analysis 

In this section, these definitions and basic facts can be found in Darchin's [6] mini-course. 

Let S be the Schwarz class, (ip, \) is a couple of smooth functions valued in [0,1] such that ip is 
supported in the shell C(0, |, f ) = {£ € R N \% < |£| < §}, X is supported in the ball B(0, f ) = {£ € 
R^IICKfjand 

OC 

For / E S' (denote the set of temperate distributes which is the dual one of S), we can define the 
nonhomogeneous dyadic blocks as follows: 

A^f := X (D)f = h*f with h = T- x X , 

A q f :=<p(2-*D)f = 2« N J h(2«y)f(x-y)dy with h = T~ x ^ if g > 0. 

where *, T~ x represent the convolution operator and the inverse Fourier transform respectively. The 
nonhomogeneous Littlewood-Paley decomposition is 

Define the low frequency cut-off by 

S q f:= \f- 

Of course, Sof = A_i/. The above Littlewood-Paley decomposition is almost orthogonal in L 2 . 
Proposition 2.1. For any f,g € S' , the following properties hold: 

A p A q f = Q if \p-q\>2, 
A,(Vi/A pS )eO if \p-q\>5. 
Besov space can be characterized in virtue of the Littlewood-Paley decomposition. 
Definition 2.1. Let 1 ^ p ^ oo and stl. For 1 ^ r < oo, the Besov spaces r are defined by 

/€B|, r ^( ^(2« s ||A g /|| LP r)"<oo 

q>-l 

and B^ x are defined by 

f e B' sup 2"||A,/||£, < oo. 

Definition 2.2. (J.-M. Bony J%$) Let f,g be two temperate distributions. The product f ■ g has the 
Bony 's decomposition: 

f-g=T f g+T g f + R(f,g), 

where Tfg is paraproduct of g by f , 

and the remainder R(f, g) is denoted by 

R(f,9) = J2 A i^i9 with A q := A g _i + A q + A q+1 . 
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Lemma 2.1. (Bernstein) Let k G N and < R\ < i?2- There exists a constant C depending only on 
Ri, i?2 and N such that for all 1 ^ a ^ b ^ oo and f G L a , we have 

Supp Tf C B^RrX) =► sup ||fl a /|| ifc < C^A^^^U/Ui*; 

a| = fc 

Supp .F/ c CtO.fliA.-RaA) C- fc - 1 A fc ||/|| i - < sup \\d a f\\ L a < C* +1 A fc ||/||i.. 

a|=fe 

Here, T represents the Fourier transform. 

A result of compactness in Besov space is: 

Proposition 2.2. Let 1 ^ p, r ^ oo, s G R and e > 0. For a// € C^°, i/ie map f t—> (f)f is compact 
from B s + r e to B s p r . 

Finally, we state a result of continuity for the composition to end this section. 

Proposition 2.3. Let 1 ^p,r ^ oo, / be open interval o/R. Let s > and n be the smallest integer 
such that n > s. Let F : I -> M safe/?/ F(0) = and F' € W n <°°(I;R). Assume that v € £* r tofces 
values in J CC /. Tften F(t>) G Bp r and there exists a constant C depending only on s, /, J and N such 
that 

iin«)iiB5„<c(i+ii«ii i =o)»iiF'ii w »,» (i) ii«i| fl . ir . 



3 Reformulation of the Original System 

In this section, we are going to reformulate (1.1)-(1.3) in order to obtain the effective a-priori estimates 
by spectral localization method. For the isentropic case (7 > 1), let the sound speed 



and denote the sound speed at a background density n by ip — ip(n). Similar to that in Ref. [IB], we 
define 

m = (^ijj(n) — %j> 

Then the system (1.1) can be reduced to the following system for C 1 solutions: 

m t + i/idivu = — u • Vm — ^-S-mdivu, 

u t + ^Vra + a = -u • Vu - Y™Vm + e, (3.1) 
e t = -VA _1 V • {/i(m)u + nu}, 

where e := V$, n(rn) = {(yl7) - ^ (^im+i/>)} 7-1 — n, is a smooth function on the domain {m|^-!-m+i/j > 
0} satisfying h(0) = 0, the non-local term VA -1 V • / is the product of Riesz transforms of /. The initial 
data (1.3) becomes into 

(m, u, e)| t=0 = (mo, Uo, e ) (3.2) 

with 

2 



m = — (jp(n Q ) - tpj , e = VA 1 (n - n). 



7' 

Remark 3.1. For any T > 0, (n,u,e) G C 1 ([0,T\ x R N ) is a solution of the system (1.1)-(1.2) with 
n > 0, then (rn,u,e) G (^([O.T] x R N ) solves the system (3.1)-(3.2) with ^ra + i/i > 0; Conversely, if 
(m,u,e) G C^^TjxR^) solves the system (3.1)-(3.2) with ^-m+ip > 0, then (n,u,e) G C 1 ([0,T]xR N ) 



2 

is a solution of the system (1.1)-(1.2) with n > 0, where n = ip^i^^m + V-0- 
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For the isothermal case (7 = 1), let n = \/A(\ nn— Inn), e = V$, then the system (1.1) is transformed 
into the following one for C 1 solutions: 



n t + \/^4divu = — u • Vn, 

u t + VAVn + ^ = u- Vu + e, (3.3) 
e t = -VA _1 V ■ {h(h)u + nu}, 

where h(n) = n(exp(A~?n) — 1) is a smooth function on the domain {n\ — 00 < h < +00} satisfying 
h(0) = 0. The initial data (1.3) turns into 

(n,u,e)| t= o = (vC4(lnn - lnn),u ,e ). (3.4) 

Remark 3.2. For any T > 0, if (n, u, e) € C 1 ([0,T]xR N ) is a solution of the system (1.1)-(1.2) with n > 0, 
then (n,u,e) e ^([O.T] x K w ) solves the system (3.3)-(3.4); Conversely, if (n,u,e) e C^tO.T] x K w ) 
solves the system (3.3)-(3.4), then (n, u, e) e (^([O,! 1 ] x l w ) is a solution of the system (1.1)-(1.2) with 
n > 0, where n = nexp(yl~2n). 

In what follows, we shall only study the system (3.1)-(3.2) and prove the main results in this paper, 
since (3.3)-(3.4) can be discussed through a similar process. 

4 Local Existence 

In this section, we first give some estimates of commutators in Besov space B^. Then, using the 
regularized means and compactness argument, we complete the proof of proposition 4.1. 

By using the first and third equation of Eq.(l.l), we get the following formulations under the variable 
transform immediately. 

Lemma 4.1. 

A ,. divet + div(ft(m)«) 

dive = him), divu= . 

n 

□ 

Applying the operator A q to (3.1) yields 

d t A q m + (u • V)A q m = — ^A^divu + [u, A q ] ■ Vm — A ? (mdivu), 
d t A q u + (u • V)A 9 u + iA,u = -ijA q (ym) + [u, A,] • Vu 

_2_iA 9 (mVm) + A 9 e, 1 ' ' 

d t A q e = -VA" 1 V • A q {h(m)u + nu}, 

where the commutator [/, g) = fg — gf. 

Multiplying the first equation of Eq.(4.1) by A q m, the second one by A 9 u and adding the resulting 
equations together, then integrating it over R N , we obtain 

lj t (\\ A « m W 2 L* + W A i u Wh) + ^\\ A «"\\h 
= \j divu(| A q m\ 2 + \A q u\ 2 ) + J ([u, A,] • VmA ? m + [u, A q ] • VuA,u) 

--^-y— J (A 9 (mdivu)A 9 m + A,(mVm)A,uJ + J A q eA q u. (4.2) 

The electric field term can be estimated as 



/ 



A q eA q u 



- / A 9 $A 9 divu 
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4 J A g $A g (drve t + div(/i(m)u) 

~ J A q eA q (h(m)u). (4.3) 



1 d ii a n9 



Note that the bi-linear spectral localization term, we have 

/ { A 3 (mdivu)A,m + A 3 (mVm)A g u} 

= - J A g m(Vm • A g u) + J [A q , m] Vm • A q u + J [A q , m]divuA g m. (4.4) 

Here, we give a lemma to estimate these commutators in (4.2) and (4.4). 
Lemma 4.2. Let s > and 1 < p < oo, the following inequalities are true: 

2^\\[f,A q ]Ag\\ LP 

{ Cc||/||B ; Jff|| B;i , /.« £ ^ s = l + 7V/p; 
^ i Cc q \\f\\ B} J\g\\ B;+ ^ f€B; A , g E B s p +\ s = N/p ; (4.5) 

[ C^II/II^HsIIb^, /eB^ 1 , g e B; A , s = N/p . 
In particular, if f = g, then 

2« s ||[/, A q ]Ag\\ LP < CcHV/llL-llffUfl;^, a > 0. (4.6) 

where the operator A = div or V, C is a harmless constant and {c q } denotes a sequence such that 

\\{c q }\\n < I- 

Remark 4.1. Similar estimates of commutators have been obtained by Chemin[3], Danchin[H 02 [6] and 
IftimiefTTj ei a/., therefore, we omit the details. 

Now, we give a proposition on the local existence and uniqueness of classical solutions to (3.1)-(3.2). 

Proposition 4.1. Suppose that (mo, Uo, eo) G B^ x satisfying ^^-mo+'ip > 0, then there exist a timeTo > 
and a unique solution (m, u, e) of (3.1)-(3.2) such that (m, u, e) G C 1 ([0,To] x W N ) with 2^— m + i/i > 
for all t G [0, T ] and (m, u, e) G C([0, T ], B° x ) n ^([0, T ], SJi 1 )- 

Proof. (Existence) Let /To = (mo,uo,eo) T G -BJ i- 

Claim 1: There exist two positive constants rf,rf' such that 

< rf < ^y^m + < »/'. (4.7) 

In fact, by the imbedding B^ 1 Co(continuous bounded functions which decay to zero at infinity), we 
know mo G Co- For any e > 0, 3Af = M(e) > such that |mo| ^ e for all |x| > M. We may choose 
e = then < |- ^ + 1/) ^ On the other hand, when a; G [-M, M], we have 

7 — 1 . -7— 1 - 7 — 1 - 

< — - — min mn + w ^ — - — mr> + w ^ — - — max mo + w. 

2 xe[-Af,Af] 2 2 k6[-m,m] 

So we obtain (4.7) only by choosing 

, f ib 7- 1 . T -| .. . r 3r5 7- 1 -1 

?7 = max < — , mm mo + tp > , 77 = mm <^ — , max mo + u> > . 

12 2 xe[-M,M] i 12 2 xe[-m,m] > 

From (4.7), we can obtain < §' ^ no ^ $" for two positive constants Note that 

'7-1 



n = no — n 
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is a smooth function on {mol^rp-mo + if) > 0}, by Proposition 2.3, we have 

IK - "-llsf,! = II'i(wo)||b|' i1 < cIKoIIb-j- 

There exists a sequence {U k } := {(n§, u§, Gq) t } such that {(rig - n,u§,eo)} G iT 2 ^ > cr, £ G Z) 
converges to (no — n, Uo,eo) in -BJ 1 and satisfies < i§' ^ n§ $J + i??'. Furthermore, we also have 

11^0 II -821 ^ ll^olls?! + 1(Uq := (mo, ug, e§) T ). We define a sequence {?7 fc } = {(n fc , u fe , e k ) T }, which 
solves the following equations: 

n\ + V • (n fc u fc ) = 

u* 1 + (u k • V)u fe + -^Vp(n k ) = e fc - ^ (4.8) 
e, fc = -VA^V • (n^u fe ) 

with the initial data 

(n k ,u k ,e k )\ t =o = (n k ,u k ,e k ). (4.9) 

It is easy to see (4.8) is a strict hyperbolic symmetric system on G = {U k \n k > 0} in the sense of 
Friedrichs. Using Kato's classical result in Ref. [13] or [15] and the L 2 - boundcdness of Riesz transfor- 
mation, we can get the following local existence result: there exist a time T k > and a solution U k to 
(4.8)-(4.9) such that 

U k G (^([0, T fe ] x R N ) with < $' k ^ n k < < for all i G [0, T k ] 

(i?j. , ^J! are two positive constants depending on k) 

and 

(,/ - n, u fe , e fc ) G C([0, T k ], H e ) n ^([0, T fc ], ff'" 1 ). 
From Remark 3.1, the sequence {?/ } = {(m fe , u fe , e fe ) T } solves the following equations: 

to^ + i/'divu' 8 = u k ■ Vm' — ^^TO^divu^ 

u t fc + $Vm k + i£ = - u fe • Vu fe - tlm'Vm' + e fc (4.10) 
= - VA- 1 V • {/i(m fc )u fc + mi fc } 

with the initial data 

(m k ,u k ,e%= = (-l_(v,(n§)-^),ug,eg) (4.11) 

satisfying 

c/ fc g ^([0, T fc ] x r n ) with < n fc < ^y- mfe + $ < Vk for a11 * e [°> 

(?7fc , ?7fc are two positive constants depending on k) 

and 

(7 fc G C([0, T k ],H l ) n C x ([0, Tfc], ff'- 1 ). 

Let [0, T£) be the maximal existence interval of above solutions to (4.10)-(4.11). Similar to the discussion 
in Ref. [15j . we have the blow-up criterion: 

T^ < 00 O limsup(||f/ t fe || L oo + ||VC/ fc || L oo) = +00 or 

for any compact subset KddG, U k (x,t) escapes K as t—*T%, 
where G = \u k \^-^-m k + ^ > oj. 

Claim 2: For t G [0, min{T*, T }), it yields \\U k {t)\\ BS>1 «S 4A , where T = A = ||Db|| B?il + 1 and 

C is a positive constant (independent of k) given in (4.18). 
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The proof of Claim 2 relics on the standard continuity method. Let 

G = \T e [0,min{T fe *,T }) : sup £(t) < 4A |, (4.12) 

where 

£(t) = \\U k (t)\\ Bli . 

In fact, we have already known 

lltfolk*,, <A <4A . 

Then from the continuity of £ (t) on [0, min{T^, To}), we may see that the set G is nonempty and relatively 
close in [0, minlT,*, To}). To show that it is also relatively open in [0, min{Tjf, To}), and hence the entire 
interval [0, min{T£, To}), it suffices to the weaker bound in (4.12) implies 

\\U k (t)\\ BSi < 2A , te [0,T] C [0,min{T*,T o }). (4.13) 

To do this, we need to make the best use of flow map. 

The flow map X(t; i, x) of u fc starting from x € l w at time i <G [0, T] can be defined as 

^ =u k {t,X(t-i,x)), X(t;t,x)\ t= t=x, 



then we have 

d -m k {t,X(t-t,x)) = (d t + u k ■V)m k (t,X(t;t,x)). 



dt 

Together with the first equation of Eq.(4.10), it is easy to get (for any t € [0,T]) 



7-1 k 



m k (t,x) + tp 

= (^m fc (0,X(0;t, a ;))+^)cxp(-l y iy t divu fe ( ? ,X( ? ;t, a ;))*). (4.14) 
Thus, there exist two positive constants 771,772 (independent of k) such that 

0< W ^m l (i,i)+^^ (t,x) e [0,T] x R N . (4.15) 

Here, we need not consider the effect of relaxation term. Therefore, by (4.2)-(4.4), Holder's inequality 
and Lemma 4.2 ( take p = 2 and s — a), we get 



|A 9 C/ fc ||| 2 



< C{||Vu fe || L oc2 2 ^(||AX||| 2 + ||A 9 u fe ||i0+c 9 2«H|u fe || s? J|m fc || B? J|AXll^ 
+ Cg 2«H|Vu fe || L oo||u fc || B? J|A 9 u fe || L2 +2 2 «n|Vm fc || L oo||AXll^l|AXllL 2 
+c,2HVm fc |Uoo||mi Bf JAXlU'+^^^^ 

+2 2 ^||A 9 (/^(m fc )u' £ )|| L2 ||A 9 e' £ || L2 }, t € [0,T\. (4.16) 

(||A 9 C/ fe ||| 2 := ||A q m fc ||| 2 + ||A 9 u fe |l! 2 + l|A 9 e fc ||| 2 , 
\\A q U k \\ L 2 « ||A 9 m fe || i2 + ||A g u fe || i2 + ||A g e fc || L2 ) 

Dividing (4.16) by (2 2 « ff || A,{/ fc || 2 2 + e) 1 / 2 ( e > a small quantity), we have 

±(2^\\A q U k \\l 2 +ey 
^ C{||Vu fc || L =o2^(||AXlli 2 + l|A g u fc || L2 ) + C9 ||u fe || B . i ||m' £ ||Bj il 
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+ Cg ||Vu' £ || LO o||u fe || s . i + ||Vm fc || L o 2^||AXll^+c 9 ||V7n fe || L oe 

x||TO fe || B?i + C9 ||m fc || B?i ||u fe || B?i +2^\\A q (h(m k )u k )\\ L 2}, te [Q,T]. (4.17) 

Integrating (4.17) on the variable t, then taking e — * and using Proposition 2.3, we obtain the a-priori 
estimate of U k : 



\u k (t)\\B Sl < l!f/ fe IU ?1 +c / ||c/ fc (0lll ?1 * 

Jo 

< {\\Uo\\B 5tl +l) + C f 11^)111^*, te[0,T]. (4.18) 



Furthermore, we have 

/•* 

sup \\u ^JUBgi UI^0||Sf ± T J-J T ^ ; sup ' ,I|J 

Set 



sup ||C7 fc (?)l|fl 2 % <Wo\\b Si + 1) + C / sup H^OIIbj,*, *e[0,T]. (4.19) 

O^c^t ' " Jo 0^?'^? 



AiW^dlObllfl^+lJ+C/ sup ||tr fc (0||i ?1 *. 

Jo o<^'<c 



Then, we have 

d 



df \i^C\l Ai(0) = ||E7b||B ?il +l, (6 [0,T]. (4.20) 



Let X(t) solves Riccati equation: 

^A-(7A 2 , A(0) = ||[/o||b 2V +1. (4.21) 

The time Tq = l/(2CAo) is less than the blow-up time for (4.21). Then by solving the differential 
inequality (4.20), we have Ai(f) < A(t) for t G [0,T ]. Solving (4.21) yields 



Therefore, we get 

||^ fc (i)||B ?il < 2A , t G [0,T] C [0,min{T,*,T }), 

which completes the proof of Claim 2. 
Furthermore, using Eq.(4.10), we have 

\\U k {t)\\ B °-i < A , t G [0,min{T,*,T o }), (4.23) 

where A is a positive constant only depending on the initial data Uq. There exist two positive constants 
fj, f] (independent of k) such that 

< fj sS ^-—?-m k (t,x) < fj, (t,x) G [0,mm{T£,T o }) xl 1 (4.24) 

From Claim 2, (4.23) and (4.24), the blow-up criterion implies < To < 7?, moreover, we have < To ^ 
MkTk - 

That is, we find a positive time To (only depending on the initial data Uo) such that the approxima- 
tive solution sequence {U k } to (4.10)-(4.H) is uniformly bounded in C([0, Tq], B$ a ) n C 1 ([0, T ], B^ 1 ). 
Moreover, it weak*-converges (up to a subsequence) to some U in L°°([0, To], B% t ) in terms of the Banach- 
Alaoglu Theorem (see Ref. [17] Remark 2 on p. 180 in Triebel, 1983). Since {U k } is also uniformly bounded 
in C([0,T ], J B^7 1 )(it weak*-converges to U t in L°°([0, T ], B a 2 7 1 )), then {f/ fe } is uniformly bounded in 
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Lip([0, T ], -BJi 1 ), hence uniformly equicontinuous on [0, T ] with the norm in B^ 1 . By Proposition 2.2, 
Ascoli-Arzela theorem and Cantor diagonal process, we deduce that 

<j)U k ->cf)U in C([0,T ],-B2,I 1 ) as k -» oo, for any 0G C c °°. 

The properties of strong convergence enable us to pass to the limit in (4.10)-(4.11). Indeed, U is a 
solution to (3.1)-(3.2). Now, what remains is to check U has the required regularity First, from U G 
C([0, To], B% 7 ) and an interpolation argument, we know U G C([0, T ], l) fo r an y (T ' < a - Furthermore, 
we have S q U € C([0, T ], Bfi) fo r an y 9 £ N. Then, combining with (4.18) (throw off the superscript 
k), we derive that S q U converges uniformly to U on [0, To] with the norm in B^ 1 . This achieves U G 
C([0,T ], J B 2 T 1 ). And then, using Eq.(3.1), we see that U t e C([0, T ], B^ 1 ), so U(t,'a:) G (^([O.To] x E^). 
By virtue of the flow map, we get < 8' ^ :L ^-m{t,x) + tp ^ 8" for (t, a;) G [0, T ] x 1^ according to 
Claim 1 (6', 6" are two positive constants). 

(Uniqueness) Let rh = mi — m 2 , u = Ui — U2, e = ei — e2 where U\ — (mi,ui,ei) and 
U 2 = {m 2l \i 2 ,e 2 ) T are two solutions to the system (3.1)-(3.2) with the same initial data respectively. 
Then U — (to, u, e) T satisfies the following equations: 

rh t + V>divu = — Ui • Vm — uVto 2 — : ^midivu — ^^-mdivt^, 

u t + tpVm + I = -ui • Vu - uVu 2 - ^miVm - ^toVto 2 + e, (4.25) 
e t = -VA^V • {H(m)ui + (h{m 2 ) + n)u}, 

where the smooth function H(ra) := h(m\) — h(m 2 ) satisfies H(0) = 0. 

Similar to the derivation of (4.18), from Lemma 4.2 (take p = 2 and s = a — 1), we obtain the following 
estimate: 

Mllijj- *S C f 11^)11^(11^1 {0\\ Bll + \\U 2 (,)\\ Bll )d, , for t e [0,T ]. (4.26) 

J 

By Gronwall's inequality, we have U = 0. □ 

5 Global Existence and Exponential Stability 

In this section, we first state a proposition on the global existence and exponential stability of classical 
solutions to (3.1)-(3.2). 

Proposition 5.1. Suppose that Uq G B 21 . There exists a positive constant 5 2 < min{^-, ^-} depending 
only on A, 7, t and n such that if 

\\U(;0)\\b Si1 <6 2 , 
then there exists a unique global solution U to (3.1)-(3.2) satisfying 

U G C([0, 00), B^) n ^([0, 00), B^ 1 ) 

and 

\M;t)\\B Sil + < Ci||^(-,0)|| Bfil exp(- Ml i), t > 0, 

where S\, C\ are three positive constants given by Proposition 5.2 and the positive constant C 2 is 
given in (5.19), U = (m, u, e) T and Ut = (m t , Ut, e t ) T . 

The proof of above proposition mainly depends on a crucial a-priori estimate (Proposition 5.2). To 
do this, we need some lemmas. 

Lemma 5.1. If{m,u,e) G C([0, T], B% A ) nC^O, T], B^J 1 ) is a solution of Eq. (3.1) for any given T > 0, 
then 
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sc 2|H| L oo(||A 9 Vm|| i 2||A ? m t || i 2 + \\A q Vu\\ L 2\\A qUt \\ L 2) + ||Vu|| L ~(||A,m t ||£ 2 

+||A,«t||i a ) + (7 - l)||m t || L ^(||A 9 divn|| L 2||A 9 mt|| L 2 + || A g Vm|| L 2 1| A g ut\\ L i) 

+(7- l)||Vm||Loo||A,mt|| L 2||A,tit|| L 2 + (||[ut,A g ] • Vm|| L 2 + A q ] ■ VmtlM 

x||A 9 m t || L 2 + (||[ut,Ag] • Vu|| i2 + ||[«,AJ • \7u t \\ L 2)\\A qUt \\ L 2 + (7 - 1) 

x(||[m t , A 9 ]divu|| L 2 + || [m, A q ]divut \\ L 2)\\ A q m t + (7 - l)(ll i m t, A 9 ]Vm|| L 2 

2 / 

+ \\[m, A q ]Vm t \\ L 2)\\A q ut\\ L 2 + -[\\A q {U{m)m t u)\\ L 2 

+ {A 1 y 1 2n^ L \\A q (m t u)\\ L 2 + \\A q (h(m)u t )\\ L 2^\\A q e t \\ L 2, (5.1) 

where H(m) = (A-f)^^ {(Aj)^^ (^m + i/i)}^ — (A-f)^^h^ 1 is a smooth function on {m\^-^-m + %l) > 
0} satisfying Tt(0) = 0. 

Proof. By differentiating the first two equations of Eq.(4.1) with respect to variable t once, integrating 
them over 1* after multiplying A q rrit, A q u t respectively, similar to the derivation of (4.16), through 
tedious but straightforward calculations, we can obtain (5.1). □ 

In addition, we give some auxiliary estimates, which are divided into high- and low- frequency cases. 

Lemma 5.2. If (m, u, e) e C([0, T], B a 2 A ) nC 1 ([0, T], B^J 1 ) is a solution of Eq. (3.1) for any given T > 0, 
then 

l|A 9 m t ||| 2 

(c^2 q \\A q u\\ L 2 + || u || L oo||A 9 Vm|| L 2 + X^|| m || ioo || A g divu|| L 2 

+ ||[u,A ? ]Vm|| L 2 + l^||[m,A 9 ]divu|| i 2)||A g m t || L 2 (q > -1); (5.2) 

2^||A,m||£ a 

-||A 9 m||l2 + ||A 9 n t || L 2 + ||ti||_L°° || A g V«|| L 2 

+ ||[u,A 9 ]-Vn|| L 2 + l_|| m ||^ ||A 9 Vm||i2 

+ I y^ll [m, A q ] ■ Vm\\ L 2 + \\A q (h(m)m)\\ L 2^ \\A q m\\ L 2 (q > 0); (5.3) 

{A^-ifi^WA-nnfv 

— ||A_in|| L 2 + ||A_in t || i 2 + ||«||l°° II A_iVtt|| i2 

+ ||[u,A_i] • Vu|| L 2 + ^— !-||m|| L oo||A_iVTO|| L 2 

+ I y^llKA_ 1 ] • Vm|| ia + ||A_i(ft(m)m)|U 2 )||A_im|| ia ; (5.4) 
2 9 l|A g e||i2 

< C{(i47)"'n^l|A,m|l^ + l|A,(Mm)m)|| L 2}||A g e|| L2 (g > 0); (5.5) 
l|A-ie|| 2 L 2 

^ (||A_ 1% || i 2 + -||A -iu\\ L 2 + CV>||A_im|| L 2 + ||A_i(tt- Vu)|| L 2 

+2zJi||A_ 1 (mVm)|| ia ) ||A_ ie |U a ; (5.6) 
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l|A g e t ||| 2 sC C(\\A g (h(m)u)\\ L 2 + n\\A q u\\ L 2)\\A q e t \\ L 2 {q ^ -1), (5.7) 

where h(m) = h'(qm)d<; - (A^^in 1 ^ is a smooth function on {ml^^-^m + tp > 0, <r G [0,1]}, 
h(Q) = and C is a uniform constant independent of A,j,t and fi. 

Proof (1) Using the first equation of Eq.(3.1), we have 

m t = — ^divu + u • Vm + — mdivu) . (5.8) 

By applying the operator A q (q ^ —1) to (5.8), integrating it over R N after multiplying A q m t , we can 
arrive at (5.2) with the aid of Holder's inequality. 
(2) Using the second equation of Eq.(3.1), we get 

ViVm = -(u t + - + u • Vu+ ^— !-mVm - e) . (5.9) 

By Lemma 4.1, we have dive = h{m). Integration by parts gives 

J A,e • VA,m = — J A q diveA q m = — J A q h(m)A q m 

= -J A q (h(m) - h(0))A q m 

= -{A^-^n^WAqmW 2 ^- j A q (h(m)m)A q m, 



where h(m) = h'(<;m)d<; — (A'y) JiiV is a smooth function on {ml^-^sm + ip > 0, <r £ [0,1]} 
satisfying h(0) — 0. By applying the operator A q to (5.9), integrating it over R N after multiplying 



AqVm, we can obtain 



^||A 9 Vm|| 2 L2 + (A 7 )-^^||A 9 m|l! 2 
< ^-||A 9 u|| L 2 + || A,u t 11^,2 + ||u||l«||A,Vu|| L 2 + ||[u, A 9 ] • Vu|| L 2 

+ 2L_l|| m || iM ||A 9 Vm|| L 2 + 2zJl||[ m ,A,] • Vm|| i 2)||VA 9 m|| L 2 
+ ||A 9 (h(m)m)|| i 2||A g m|| L 2. (5.10) 
By Lemma 2.1, we know 

||A,Vm|| L 2 «29||A 9 m|| L 2 (q > 0), 

so (5.3) follows. For the low frequency (q = —1), by (5.10), we get (5.4) immediately. 

(3) By applying the operator A q to both sides of dive = h(m) (q ^ 0), integrating it over after 
multiplying A 9 dive, we reach (5.5) in virtue of Holder's inequality For q = —1, using the second equation 
of Eq.(3.1), we can obtain (5.6). 

(4) From the L 2 - boundcdness of Riesz transform, we can derive (5.7) directly from the last equation 
ofEq.(3.1). □ 

The crucial a-priori estimate is comprised in the following proposition. 

Proposition 5.2. There exist three positive constants 81, C\ and [i\ depending only on A,j,n and r 
such that for any T > 0, if 

sup (||I/(-,t)|| Bfl + ||I/ t (-,t)|| B j T i)<5i, (5.11) 

then 

\\U(;t)\\B Sl + \\U t (;t)\\B'-^Ci\\U(;0)\\B Sl eM-^t), t€[0,T\. 
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Proof. From the a-priori assumption (5.11), we have 



sup 

0<t<T 



(\\u(;t)\\ wl ,oo + \\u t (;t)\\ L ~) <c«i. (5.12) 



To ensure the smoothness of functions h(m),H(m) and h(m), we may choose 
< S 1 sC (7 ^ 1)c , then 



1— !^m(i,a;) + V > I > 0, (i,x)€[0,T]xl N 

and 

l^l<;m(t,x)+$ > | > 0, ?e [0,1], (t.a;) G [0,T] x R N . 
From (4.2)-(4.4), we set 

Ji,, := ||Vu||£oo(||A,m||i 2 + ||A,u||i 2 ) + 2||[u,A g ]-Vm||£2||A g m|U» 
+2||[u,A g ] • Vu|| L 2||A 9 u|| i2 + ( 7 - l)||Vm|| L »||A,m|| L 2||A g u|| L 2 
+ ( 7 - 1)|| [m, A 9 ]Vm|| i 2||A 9 u|| L 2 + (7 - 1)|| [m, A g ]divu|| L 2 1 |A g m|| L 2 

+ |||A 9 (^(m)u)|| i 2||A 9 e|| L 2. (5.13) 

l2, q denotes the right-hand side of inequality (5.1). For the proof of Proposition 5.2, we are going to 
divided it into the following two lemmas. 

Lemma 5.3. (q ^ 0) There exists a positive constant \i 2 depending only on A,j,t and n such that the 
following estimate holds: 

2^-i)|{ 2 2 «(||A 9 m||| 2 + \\A q u\\l 2 + |||A,e||i a ) + (||A g m t ||| 2 

+ 1^^1112 + lllA^Hi,)}' + ^ 2 (2^\\A q U\\ L2 + 2^-V\\A q U t \\ L2 ) 
< C{2^(\\U\\ W ^ + \\U t \\ L ~)(20\\A q U\\ L 2 + \\A q U t \\ L 2) + c q (\\U\\ Bli 
+11^11^-011^1^ +2^\\A q (h(m)u)\\ L2 +2^\\\A q (n(m)m t u)\\ L , 
+ \\A q (m t u)\\ L 2 + \\A q (h(m)ut)\\ L 2 + \\A q (h(m)m)\\ L i)}, (5.14) 

where C > is a harmless constant depending only on A, 7, r and h. 

Proof of Lemma 5.3. Combining (4.2)-(4.4), Lemma 5.1 and 5.2, we have (q ^ 0) 



|{2 2 «(||A 9 m|||2 + ||A g u||| 2 + ±||A,e||£ a ) + (||A g m t ||| 2 + ||A,u t || 
+ ±||A,e,||£ 2 ) } + M2 2 i\\A q m\\ 2 L2 + ^\\A q u\\ 2 L2 + ^\\A q e 



2 



+/33||A 9 m t ||| 2 + ^||A 9 u t |||2 + |||A 9 e t |||2 



^ 2 2 n hq + J 2> , +/3 1 C2«(l||A ? u|| L 2 + ||A g u t || L2 + ||u|Uoo||A g Vu|| L 2 



+ || [u, A 9 ]Vu|| L 2 + X_I|| m || L ^||A 9 Vm|| L 2 + l_I||[ m ,AjVm||i2 

+ \\A q {h{m)m)\\ L .) || A g m|U 2 + ^2«((A 7 )-5n 2 i 2 1| A q m\\ L . 

+ || A q (h(m)m)\\ L 2^ ||A 9 e|| L2 + /? 3 (cWl|A g u|| L 2 + ||u|| L » ||A g Vm|| L 2 

+^=-||m|| L oc||A,divu|| 

L 2 H — II [m, A 9 ]divu|| L 2 + || [u, A g ]Vm|| L 2 
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P±-(\\A q h(m)u\\ L 2 +n||A,u|| L2 )||A,et|| L2) (5.15) 
where these positive constants (3i,p2, Pz and (5± satisfy 

respectively. We introduce them in order to eliminate quadratic terms in the right-hand side of (5.15). 
There are no quadratic terms in Ji i9 and I 2 , q . By Young's inequality, the first quadratic term can be 
estimated in this way: 

/3 1 (7i2^||A 9 u|| L2 ||A 9 m|| L2 

T 

V If} T v 

< ^2 2 ^||A 9 u||i 2 + \M2^\\A q m\\l 2 

< ^IIA^ulli. + ^x^HAHIi- 

The remainder quadratic terms in the right-hand side of (5.15) are estimated similarly as follows: 

/JiC2«||A,u t || L2 ||A,m|| i2 

< ^!||A 9 u t ||| 2 + i/3^2 2 «||A 9 m||| 2 

< illA^n^ + ^^iiAHl^; 

^(A 7 )-5n^2«||A 9 m|| L2 ||A 9 e|| i2 

< ^^IIA.HIi. + ^IIA.elli, 



i/3 1 ^||A g m||| 2 + ^2 2 «||A gt 



C , ^/3 3 2 9 ||A 9 u|| L2 ||A 9 m t || i2 
< ^!/332 2 «||A 9 u|| 2 2 + i/3 3 ||A g mJ| 2 
= ^2 2 «||A 9 u||i 2 + i/3 3 ||A 9 m t || 2 L2 ; 



,S4C1|A 9 u|| i2 ||A 9 e t || i2 
< !^!/J 4 2 2 «||A g u||| 2 + ^/3 4 ||A 9 e t ||| 2 
= l2 2 «||A 9 u||| 2 + ^/3 4 ||A 9 e t || 2 2 . 



Then (5.15) becomes into 



|{2 2 «(||A 9 m||| 2 + ||A,u||£ a + i||A 9 e|| 2 2 ) + (||A 9 m t || 2 2 + ||A,u t | 
+ 



2 

L- 



A 9 e t ||| 2 ) } + ^2 2 «|| A 9 m|| 2 L2 + ^2 2 «|| A 9 u||i 2 + ^2 2 «|| A 9 e||i 2 
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+ ^\\A q m t \\l 2 + ±||A,u t ||£ a + §l|A 9 e t ||| 2 
< 2^h. q + I 2 . q + /3 l( 72«(||u|| L oo|| A 9 Vu|| i2 + 2_!|| m || ioo || A 9 Vm|| i2 



+ || [u, A 9 ]Vu|| i2 + l_||[ m ,A 9 ]Vm|| L 2 + \\A q {h(m)m)\\ L 2^\\A q m\\ L 2 
+ ^29||A g (h(m)m)|| L 2||A,e|| L 2+/3 3 f||u|| LO c||A g Vm|| L 2 + li: -\\m\\ L o 



x||A 9 divu|| L 2 + || [u, A 9 ]Vm|| L 2 + 2_||[ m , A 9 ]divu|| i 2 J || A q m t \\L 2 
(3iC 



n 



A 9 /i(m)u|| L 2||A 9 et||z,2. (5.16) 



Dividing (5.16) by 

{2 2 ^(||A g m|| 2 L 2 + ||A g uf L 2 + ±||A,e||£ a ) + (||A,m t ||£ 2 + ||A g u 4 || 2 L2 + ^||A g e t ||| 2 ) }* 

and multiplying (5.16) by the factor 2 gt -' J ~ 1 \ we get (5.14) immediately with the help of Lemma 4.2, 
which completes the proof of Lemma 5.3. □ 

For the case of low frequency (q = — 1), we also have the following a-priori estimate in a similar way. 

Lemma 5.4. (q = —1) There exists a positive constant ^3 depending only on A,j,t and n such that the 
following estimate holds: 

T |{I(||A_ 1 m||i a + l|A-m||£ 2 + i||A_ ie ||£ a ) + (\\A- imt \ 



< C 



1 d J 1 f 1 ! ||2 \ \ 1 ' ' 1 ■ 1 1 \ ' \ 1 / 1 1 \ 

A_ lU J 2 2 + i||A_ ie J 2 2 )} 2 + M2 (1.||A_ 1 C7|| L 2 + _2 T ||A_ 1 t7 t || L 2^ 

{2^1(11^^.- + ll^tlkoc) (^IIA-^H^ + HA.^.iu.) + c- 1 (\\u\\ Bll 

+ 11^11^-011^11^ + ^||A_!(/i(m)«)|| La + ^ T (||A_ 1 (W(m)m t u)|| La 
+||A_i(mt«)|| ia + ||A_i(/»(m)tit)|| La + || A_i(/i(m)m)|| L 2) }, (5.17) 

where C > is a harmless constant depending only on A,j,t and n. 

Summing (5.14) on q £ NU{0} and adding (5.17) together, according to Proposition 2.3 and a-priori 
assumption (5.11)-(5.12), we get the following differential inequality: 

|q + M4 (||C7|| b?i + \\U t \\ BVil ) < C5 1 {\\U\\ B , 1 + \\U t \\ B „-,), (5.18) 



where 



Q: = 2^- 1 ){2^(||AHl! 2 + l|A 9 u|||2 + i||A 9 e||£ 2 ) 

+ (||A 9 m t ||| 2 + ||A 9 u 4 ||£ 2 + i||A 9 e t ||| 2 )}" 



n 

and the constant 114 depends only on A, 7, r and n. Furthermore, it is easy to show that Q satisfies 

Cs(\\U\\ BSl + ||^|| B j-i) < Q < C 4 (\\u\\ BSi + \\u t \\ B ,-. 

for two positive constants C 3 , C 4 . Choosing 5i = min{^, ^^ c }, we complete the proof of Proposition 
5.2 with Hi = f. □ 
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Proof of Proposition 5.1. From the assumption 

\\U(;0)\\ Bll ^S 2 , 

by Proposition 4.1(local existence), we can determine a time T\ > (Ti < T ) such that 

M-Mbs^ + \\Ut(;t)\\ B ^ < C 2 \\U(;t)\\ BSi < 2C 2 S 2 (5.19) 

for all t € [0,Ti]. 

Claim: One can choose a positive constant (52 satisfying <5 2 < min{^-, ^-} to ensure 

II^(-^)I|b 2 V + \\Ut(;t)\\ B ^ < Si, for all t e [0,T ]. (5.20) 

Otherwise, we may assume that there exists a time T2 (Ti < T 2 < To) such that (5.20) is satisfied for all 
t S [0,T 2 ) and 

\\U(;T 2 )\\ BSi + ||^(-,r 2 )|| BJ -i - fc, (5.21) 

because (5.20) is satisfied as t G [0, T{\ for such a choice of S 2 . By Proposition 5.2, for all t G [0, T fe ] (T fc — > 
T 2 — ) 

+ II^M)!^ 1 <C 1 ||[/(-,0)|| S j i exp(- Ml i). 

In particular, 

\\u{-,T k )\\ Bli + \\u t {-,T k )\\ B ^ < c 1 ||[/(.,o)|| B?i ex P (- Ml T fc ) 

s$ C 1( 5 2 exp(- Ml T fe ) < dS 2 . 

By the continuity onte [0, T ], we get 

\M;T 2 )\\ BSi + \\U t (;T 2 )\\ BVil < < S 1 , 

which contradicts (5.21). So, (5.20) holds. From Proposition 4.1 and 5.2, using the boot-strap argument, 
we can prove Proposition 5.1. □ 

From the imbedding property in Besov space B 21 , (m, u, e) G C 1 ([0, 00) x R N ) solves (3.1)-(3.2). 
The choice of 5i is sufficient to ensure m + ?/> > 0. According to Remark 3.1, we know (n, u,e) G 
C 1 ([0,oo) x R w ) is a solution of (1.1)-(1.3) with n > 0. Furthermore, we may attain Theorem 1.2. 

Finally, we show the exponential decay of the vorticity. 



Proof of Theorem 1.3. When N = 2 and 3, the curl of the velocity equation in Eq.(l.l) gives 
d t uj 4 

Then, we may get 



d t co + —u + u • Vw — oj ■ Vu = (in particular, u ■ Vu = when N = 2). 

T 



1 d .. . ,,, 1 „ , 
- s l|A^||l 2 + -||A gW ||| 2 

< C 5 (||Vu||Loo||A g a;|| L 2 + \\u>\\ L oc || A 9 Vu|| i 2 

+c,||Vu|| B? -i||w|| B? -i)||A,w|| L2 . (5.22) 



Dividing (5.22) by ||A 9 w|| i 2 and summing it on q ^ —1 (q G Z) after multiplying the factor 2 q ^ from 
Theorem 1.2, we have 

- — lk(-,*)|| B? -i + ;||w(-,t)|| B? -i 



2d* 1 "2,1 r 

< C 6 ||ll(-, t)||B ?>1 

s$ C 5 C (\\(n - n,u ,e )\\ BSi )\\u)(-,t)\\ B „-i 
s? C 5 C min { J , 1 - } 1 1 u> ( ■ , t) \ \ B « - 1 



< ^llK-,t)|| B?7 i. (5.23) 



1 

Therefore, we obtain the exponential decay of |jw(-, i)|| B „-i. □ 
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